Enhancing Tc in ferromagnetic semiconductors 
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We theoretically investigate disorder effects on the ferromagnetic transition ( "Curie" ) temperature 
Tc in dilute IIIi-j^Mnj^V magnetic semiconductors (e.g., Gai_i,Mna,As) where a small fraction {x « 
0.01 — 0.1) of the cation atoms (e.g., Ga) are randomly replaced by the magnetic dopants (e.g. Mn), 
leading to long-range ferromagnetic ordering for T <Tc. We find that Tc is a complicated function 
of at least eight different parameters including carrier density, magnetic dopant density, and carrier 
mean free path; nominally macroscopically similar samples could have substantially different Curie 
temperatures. We provide simple physically appealing prescriptions for enhancing Tc in diluted 
magnetic semiconductors, and discuss the magnetic phase diagram in the system parameter space. 
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Diluted magnetic semiconductors (DMS), which are 
ferromagnetic for T < Tc, are of great fundamental inter- 
est because they provide a unique example of interplay 
among magnetism, disorder, carrier dynamics, and trans- 
port properties. Starting with the early report of carrier- 
induced global ferromagnetism first in Ini-^^Mn^^As [1] 
and later [2] in Gai-ajMn^jAs, with the typical magnetic 
impurity concentration x ~ 5%, there have been many 
studies reporting ferromagnetism in DMS systems as dis- 
parate as Gei-xMux [3], Gai_a;Mn2;N [4], Gai_a;Mna;P 
[5], Ti02-Co [6], Sn02-Co [7], and ZnCrTe [8] among 
others. There is a great deal of current research activity 
in the subject, and new DMS materials with novel mag- 
netic properties are likely to continue appearing in the 
near future. In spite of this enormous activity, there is 
no current consensus on the basic magnetic model un- 
derlying DMS ferromagnetism. In fact, there is strong 
disagreement in the literature about the nature of the 
carrier ground state participating in the ferromagnetism 
— in particular, whether these carriers are free carriers in 
the parent-semiconductor bands or are impurity (or de- 
fect) band carriers created within the fundamental band 
gap of the semiconductor by the magnetic dopant atoms 
is a matter of considerable controversy in the literature 
[9]. Experimentally the situation turns out to be quite 
complex as the ferromagnetic properties (in particular, 
the Curie temperature Tc) seem to depend very sensi- 
tively [10] on the materials growth and processing (e.g. 
annealing) conditions, and Tc in nominally identical sam- 
ples could differ substantially depending on the precise 
details of sample preparation and processing. 

In this letter, we develop a theory for the prediction 
of Tc in DMS systems by focusing on the most-studied 
DMS system, viz. Gai_a;Mna;As, as a function of various 
system parameters (to be discussed below). Our work, 
which includes spatial disorder of the system (by virtue of 
the random locations of the magnetic Mn atoms), leads 
naturally to the conclusion that Tc is a complicated (and 
in general unknown) function of the system parameters; 
the number of independent variables (at least eight as 
discussed below) determining Tc, even within the min- 



imal zeroth order effective Hamiltonian approach, are 
sufficiently large that nominally macroscopically identi- 
cal samples may very well have significant variations in 
Tc, as is indeed observed. Some of the system parame- 
ters (e.g. the nature of spatial disorder, various defects 
present in the real samples, the impurity scattering po- 
tential, etc.) are generally unknown as a matter of prin- 
ciple, and therefore a precise quantitative prediction for 
Tc as a function of the known system parameter (i.e. the 
Mn concentration x and the carrier density Uc), as is of- 
ten done in the literature, could be quite meaningless 
because the same values of x and nc may lead to dif- 
ferent Tc values in different samples depending on the 
details of various defects/impurities/disorder in the sys- 
tem (which can be influenced, for example, by sample 
annealing). Among the various possible defects. As an- 
tisites and Mn interstitials are known to be important. 
In addition, likely correlated clustering of Mn atoms (in- 
stead of uncorrelated random positioning at Ga sites) 
is also thought to be significant [11]. Experimentally, a 
strong dependence of Tc on the sample conductivity has 
been found with Tc typically increasing with conductiv- 
ity; samples with the highest (lowest) conductivity in- 
variably have the highest (lowest) TcS, but whether this 
dependence arises entirely from the higher conductivity 
samples having higher carrier densities (as is commonly 
assumed in the literature) or longer mean free paths 
(MFP) (or equivalently, higher mobility) or a combina- 
tion of the two is presently unknown. We therefore see 
that Tc = Tc{ni,nc, I, Jo; Nas, Nj, Nc, Jaf ■ ■■) is a func- 
tion of at least eight different parameters due to the com- 
plexity of the problem. The three obvious parameters, 
which have been discussed widely in the literature, are 
the active magnetic moment density rii due to substitu- 
tional Mn dopants, the carrier density ric, and the carrier 
MFP I (as obtained from the dc-conductivity) , and in this 
paper we mostly focus on the dependence of Tc on these 
three parameters. The other relevant parameters, e.g. 
Nas (the As antisite defect density), iV/ (the Mn intersti- 
tial defect density) , Nc (a set of parameters defining the 
clustering of Mn atoms or the correlation in their spatial 
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positions), and Jaf (the direct short-range Mn-Mn an- 
tiferromagnetic exchange interaction), are quantitatively 
important, but essentially unknown (either experimen- 
tally or theoretically), and their explicit inclusion in the 
theory for a quantitative comparison with the experimen- 
tal results is hence not particularly useful. However, we 
assume that the influence of these parameters can be in- 
cluded qualitatively in the theory by appropriately ad- 
justing the parameters rij, Tic, I, and most importantly, 
the effective ferromagnetic coupling strength Jq which 
sets the overall energy scale in the problem since we ex- 
press Tc in the units of Jq. We find that even with just 
three independent parameters {jii, Uc, I), the problem is 
quite rich leading to many subtle possibilities. 

We use the standard RKKY-Zener effective magnetic 
model [12] for describing the coupled carrier-local mo- 
ment system since we are interested in the so-called 
"metallic" DMS regime with itinerant carriers where Tc 
is maximized. (The insulating strongly localized DMS 
regime can also be ferromagnetic, but the nature of the 
insulating DMS ferromagnetism, with typically rather 
low Tc, is fundamentally different [13] from the metallic 
regime with higher Tc being considered in this Letter.) 
We note that while it is somewhat crude to characterize 
the metallic DMS regime by a MFP I extracted from the 
conductivity a = rice^T/mc (and I = vpT, where r is 
the transport relax;ation time, rric the carrier mass, and 
vp the carrier Fermi velocity), the extracted value of the 
MFP I is rather short (typically around one to a few lat- 
tice constants). Following our earlier work [14] we employ 
an effective magnetic model Hamiltonian for interacting 
impurity moments where the carrier degrees of freedom 
have been integrated out: 

H = Y, JF{Vi,)Si ■ S, + Y^JAF^i ■ S„ (1) 

where the subscripts F (AF) denote the ferromagnetic 
(antiferromagnctic) part of the impurity spin (S^) inter- 
action and the prime in the second term implies that the 
sum is restricted to nearest neighbors since the short- 
range antiferromagnctic interaction only couples nearest- 
neighbor Mn atoms (if there are any). We do not explic- 
itly include the direct antiferromagnetic exchange term in 
the results shown in this paper assuming without any loss 
of generality (within our theoretical scheme) that Jaf ef- 
fectively modifies (in fact, reduces) the carrier-mediated 
ferromagnetic interaction. Since the interaction strength 
is the basic energy parameter in our effective magnetic 
Hamiltonian, it makes little sense to keep two unknown 
coupling parameters. We write Eq. (1) with the direct 
antiferromagnetic exchange term for our general qualita- 
tive discussion of the DMS phase diagram as described 
below. 

The first term in Eq. (1), the carrier- mediated ferro- 
magnetic inter-impurity interaction, is of the RKKY form 



in our effective model: Jpij) = Jo[(2fci?r) cos(2fci?r) — 
sm{2kFr)]/ {kprY , where Jq is the fundamental ferro- 
magnetic coupling parameter in the problem (which im- 
plicitly includes all materials and band structure infor- 
mation about the system) and kp the carrier Fermi mo- 
mentum (and r = \Yij\ the spatial separation between 
randomly located substitutional Mn pairs in the GaAs 
lattice). We note one important feature of Jf (r) : for high 
carrier density, Uc > rii, the oscillatory aspects of RKKY 
interaction come into play, potentially suppressing DMS 
ferromagnetism — the details of this suppression are an 
important topic of this work. A straightforward Weiss 
mean-field treatment of the RKKY interaction, neglect- 
ing all spatial disorder effects and assuming a continuum 
virtual crystal approximation (VGA), was first carried 
out a long time ago (and has recently been rediscov- 
ered in the DMS context) leading to: T^-^^ oc riinl^^. 
This TY'^^ (with appropriate quantitative modifications 
arising from band structure effects) has been extensively 
(and uncritically in our view) used in the DMS literature 
for explaining and predicting Tc in DMS systems. In the 
current work we include in our calculation of Tc spatial 
disorder (thus relaxing the continuum VGA) and finite 
MFP, which are both significant in DMS materials. 

First we qualitatively discuss the DMS phase diagram 
as a function of the variable rii, He, and I, noting that 
the standard VGA implies that the system is a ferro- 
magnct for all values of Uc and rii with tJ^'-^"^ increasing 
monotonically with increasing impurity (rij) and carrier 
[uc) density. It is essential to include MFP effects in 
the theory. The RKKY interaction has been calculated 
in the presence of resistive scattering earlier in the lit- 
erature [15], and the modified RKKY interaction in the 
presence of a finite MFP / has the form J_F(r; I) = J Fir) 
for r < / and JF{r;l) = Jq cos(0(r))/(2A;Fr)^ for r » 
where Jpir) is the standard RKKY formula and (p{r) is a 
completely random function of r. The inclusion of short- 
range direct antiferromagnetic exchange [Eq. (1)], RKKY 
oscillations, and scattering/transport MFP effects per- 
mits us to obtain the qualitative T = DMS phase dia- 
gram as a function of three length variables Ac ~ Uc^^^; 

— 1/3 

To ^ ' ; I. The schematic phase diagram depicted in 
Fig. 1 assumes Ac, To, and / to be completely independent 
variables (which they cannot be in real systems). In each 

RKKY"SG"g FMI i FM II ^ ^ 
$5 Tj, increasing ; Tj, decreasing 

"RSO" g FMI : FMII |RKKY"SG'^ ^ 
(with Jaf) Tj, increasing i decreasing Jj; ^ 

, . disordered "RG" ^ FM I : FM , 
(c) 1 \ I 

^ T(, increasing ; T^, saturated 
FIG. 1. DMS phase diagram. Note that the dashed lines 
between FM phases are not phase transition. 
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case shown in Fig. 1 wc assume the other variables to be 
fixed at some reasonable "optimum" values (which may 
not always be experimentally possible). 

In Fig. 1(a) the system is an "RKKY spin-glass (SG)" 
system at very small values of Ac (equivalently very large 
values of ric). This arises from the frustration induced by 
the RKKY oscillations which are dominant in the high 
kp {nc ^ rii) limit (but are essentially absent in the 
dilute Tie -C rii limit). As Ac increases (i.e. ric ~ Aj^ 
decreases) the "RKKY SG" phase gives way to a fer- 
romagnetic phase. The Curie temperature of this phase 
increases with decreasing density until the optimum Tc is 
reached, after which Tc begins to decrease with decreas- 
ing carrier density. In Fig. 1(b) we depict the DMS phase 
diagram in magnetic impurity density rij (~ r^^) assum- 
ing I and Ac to be fixed at reasonable optimal values. For 
small ro (large rij) the nearest- neighbor direct antiferro- 
magnetic exchange between the Mn moments becomes 
important and competes with the carrier-induced RKKY 
interaction, leading to a random spin-ordered (RSO) non- 
ferromagnetic phase (which may or may not be a SG 
phase). For larger tq we enter the ferromagnetic phase 
with Tc (ro) similar to the profile in Fig. 1(a); even- 
tually we find a phase transition to the RKKY SG phase 
for very large rg where Uc ^ n^. Thus, except for the 
antiferromagnetic exchange induced "RSO" phase in the 
high impurity concentration limit in Fig. 1(b), the phase 
diagrams of Fig. 1(a) and (b) are essentially mirror im- 
ages of each other as one shifts from ric ^ rij {rii ^ ric) to 
Tic *C ni {uc » rii) in Fig. 1(a) [(b)] respectively. Finally, 
in Fig. 1(c) we discuss the phase diagram as a function 
of the MFP for fixed rii and ric values. For very small 
I, the carrier-induced ferromagnetic RKKY coupling is 
suppressed, and the inter-impurity interaction has ran- 
dom sign, leading to a type of (non-RKKY) disordered 
"random glassy" (RG) non-ferromagnetic ground state 
which, with increasing Z, should make a phase transition 
to the ferromagnetic phase. As one increases I further the 
ferromagnetic phase should initially be enhanced (i.e. Tc 
rises with /) with an eventual saturation of Tc determined 
by the precise values of ric and rii. Thus it is readily ev- 
ident that the 'best' technique to enhance Tc would be 
to increase the MFP sitting at fixed optimal values of ric 
and rii. Although this is perfectly reasonable as a mat- 
ter of principle, it may be difficult to increase the MFP 
without affecting ric and n^. 

Following this qualitative introduction to the DMS 
phase diagram on the basis of the simple effective Hamil- 
tonian approach, we now consider the quantitative de- 
pendence of Tc{nc, Hi, I) in DMS systems focusing on the 
well-studied Gai-^^Mn^jAs system. To do this we carry 
out a thermal lattice mean field treatment (treating ex- 
actly the spatial disorder of random Mn locations at Ga 
substitutional sites in the zinc blende GaAs lattice) of 
the effective Hamiltonian as described in [14]. Such a 
treatment avoids the physically unrealistic assumptions 




FIG. 2. Calculated Tc (solid lines) as a function of nc/rn 
(a) for different x values and fixed MFP I /a — 5, and (b) for 
different MFP l/a and fixed x = 0.07. Dashed lines indicate 
the results from VGA. 

of the continuum VGA involved in deriving Tc, and 
should be an excellent approximation for obtaining Tc 
because of the very large coordination number in the fee 
zinc blende GaAs lattice structure. As emphasized by 
us elsewhere, this lattice mean field theory (LMFT) can 
essentially be carried out for infinite size systems, thus 
avoiding the finite size complications inherent in direct 
Monte Garlo simulations which are not particularly well- 
suited for rapidly determining Tc. 

In Fig. 2 we present our theoretical results for 
Tc{nc,ni,l). As noted earlier, we have ignored the di- 
rect short-range antiferromagnetic exchange interaction 
in obtaining these results. We also ignore interstitial de- 
fects (and antisite As) in this calculation; we assume that 
Ui and ric are the effective active local moment and hole 
density respectively, which already incorporate various 
defect effects. 

In Fig. 2(a) we show the calculated Tc as a function of 
ric/rii for several values of x in a "highly" metallic system 
(l/a = 5, where a is the GaAs lattice constant), whereas 
in Fig. 2(b) we show Tc{nc/ni) for a fixed value of x 
(= 0.07) but for several different values of l/a. We show 
the simple continuum VGA result in each case for the 
sake of comparison. It is obvious that the simple theory 
of TY'~^^ is qualitatively incorrect for large ric /rii where 
Tc actually reaches a maximum and then decreases with 
increasing carrier density (due to the frustration inherent 
in the RKKY oscillations playing a role for ric/rii > 1 or 
equivalently kpro > 1) in contrast to the erroneous claim 
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(made extensively in the literature) that Tc{nc) oc ric 
(as obtained from continuum VGA) would continue in- 
creasing monotonically with carrier density. It is, how- 
ever, important to note that our results presented in 
Fig. 2 indicate a fairly large regime of (qualitative and 
even semi-quantitative) validity of the simple VGA with 
the appropriate numerical modification of Jo which is an 
adjustable parameter merely setting the scale of energy 
in our theory. In particular, the simple VGA theory re- 
mains valid up to ric/n^i 0.5, and perhaps even up to 
Tic/ rii ~ 1-0 depending on the Mn content (i.e. x). The 
optimum ric/rn value where Tc is maximum decreases 
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FIG. 3. Calculated rc(Z) for different iic/ni and fixed 
X = 0.07. Note tY^'^/Jo =32.6, 41.6, 52.2, 63.1 for n^/m =1, 
2, 4, 7, respectively. 

as Mn concentration increases, and this indicates that, 
for a given value of the effective coupHng Jq, can- 
not really be arbitrarily increased just by increasing the 
carrier density (although increasing carrier density by co- 
doping with non-magnetic impurities [16] should enhance 
Tc somewhat), but typically Tc is optimum for ric/ni ~ 1. 
This suggests that the current popular wisdom of trying 
to enhance Tc for GaMnAs (and other DMS materials) 
simply by increasing carrier density would not work much 
beyond Tc ~ 300K (since the current maximum Tc is 
around 170K). 

Our results for Tc{l) presented in Fig. 3 indicate one 
possible strategy for enhancing Tc- As can be seen in 
Fig. 3, Tc increases monotonically with increasing MFP 
eventually saturating at the maximum possible Tc for a 
given value of Uc/ni (which is somewhat below the corre- 
sponding TY^^). Thus a clear strategy to enhance Tc is 
to optimize nc/ui for a given Mn content such that one 
is at or near the optimum carrier density (i.e. near the 
maximum in Fig. 2), Uc/ni '~ 0.5 — 2 depending on x, 
and then to enhance the carrier MFP by reducing scat- 
tering effects through a systematic improvement of sam- 
ple quality. Experimentally, it is now established [17] 
that enhancing conductivity by improving sample qual- 
ity (e.g. via annealing) can substantially increase Tc, but 
the conductivity a depends both on the carrier density 
Uc and the MFP Z, and it has almost universally been 
assumed that the increase of Tc due to enhanced a arises 
entirely from the increasing carrier density, whereas we 
find that Tc improvement arises both from increasing Uc 
and /. 

Finally, in Fig. 4 we show a direct comparison between 
our theory and recent GaMnAs experimental results from 
several different groups [10,17-19]. For each set of re- 
sults in Fig. 4 we have extracted Tc{nc,ni,a{l)) from 
the relevant experimental work as described in the fig- 
ure captions. It is obvious that the experimental results 
are well-described by the theory; the agreement can be 
made essentially exact by slightly adjusting rii and/or 
by choosing slightly different Jq for different values of x, 
both of which may be reasonable since Mn interstitials 
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FIG. 4. (a) Comparison of Tc with experimental results of 
Ref. [18] and [19] (inset) as a function of x. (b) Comparison of 
Tc with experimental results of Ref. [10] as a function of hole 
density and those of Ref. [17] as a function of conductivity 
(inset). 

(whose density may very well be a function of x) are likely 
to affect the value of effective coupling by modifying the 
number of active Mn moments participating in global 
ferromagnetism and by introducing some direct Mn-Mn 
antiferromagnetic exchange coupling. It is also evident 
from Fig. 4 that the continuum VGA, while qualitatively 
reasonable in some regimes of the parameter space, does 
not provide a good quantitative description for the ex- 
perimental results. 

In conclusion, we have discussed the DMS phase dia- 
gram, taking into account effects of carrier-mediated fer- 
romagnetic and antiferromagnetic coupling between the 
impurity moments as well the frustration arising from 
the RKKY oscillations, finding that ferromagnetism is 
only one of four distinct magnetic phases possible in the 
disordered system. We have also carried out a detailed 
theoretical Tc calculation as a function of magnetic impu- 
rity concentration, carrier density, and conductivity (i.e. 
MFP) including full effects of spatial disorder and ran- 
domness, finding that the maximum Tc is obtained for 
an optimum carrier density Uc/ni ~ 0.5 — 2 depending 
on the Mn concentration. We have shown that for Tc to 
be enhanced at given values of Uc and n^, one needs to 
increase the MFP as much as possible. 
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